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I. INTRODUCTION 



Consider the two dimensional vector ordinary differential equation: 
v = A v + sgn (rv 1 + sv 2 - c) , (1) 

2 2 

where r + s > 0 , is some non-zero constant vector, and A is a 

constant matrix whose characteristic polynomial is 

A + 2yA + (y + go ), y>0 

Equations of this type commonly arise in systems involving ideal relays. 

We refer to the curve rv^ + sv^ - c = 0 as the switching line in these 
applications, and to the cases c ^ 0 and c = 0 as biased and unbiased, 
respectively. In this paper, we investigate necessary and sufficient 
conditions under which (1) will have periodic (limit cycle) solutions. 

We shall show that existence of a limit cycle in the unbiased case is a 
necessary and almost sufficient condition for existence of a limit cycle 
in the biased case, however, that not all biasings (i.e., values of 
c ^ 0) will produce cycles. We also develop an extremely powerful and 
interesting geometrical interpretation of the necessary and sufficient 
conditions, and use this interpretation to demonstrate a bifurcation 
into two cycles for the same value of c under certain circumstances. 

Equations of this type have been considered by Fleishman [2]> 

it 

Davis and Fleishman [l]> and Flugge-Lotz [3], among others. However, Davis and 
Fleishman consider primarily the special cases when only even order 

ti 

derivatives are present, and Flugge-Lotz does not give a complete 



1 



characterization of the limit cycles. Also, neither author considers the 
question of multiple limit cycles in the general case. 

For the purposes of our investigation, it is convenient to introduce 
a non-singular linear transformation and time scaling (Appendix A) that 
reduces (1) to: 

x 1 = -28 - (8 2 + l)x 2 + a 1 sgn (x 2 - c) , 

x 2 = x x + a 2 sgn (x 2 - c) , (2) 

2 2 

where 8 > 0 and a l + a 2 > ® 

This now represents a damped system, forced by a relay switching along 
the transformed switching curve = c . It can be shown that a 

necessary and sufficient condition to avoid the phenomenon called 
"chatter” in (2), and to insure that a solution will exist at each point 
in the (x^,x^) plane, is ^ 0 . We now make that restriction. 

Furthermore, (2) involves only odd functions of x^ and x^ ; hence we 
can assume, without loss of generality, that c ^ 0 
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II. THE PHASE PLANE 



Our analysis of (2) begins with the observation that, viewed 
separately in the regions > c and < c , (2) reduces to a linear, 

constant coefficient system. Elementary methods then yield the solutions 



: k * — g £ 

x^(t) = - c^e (3 cos t + sin t) + c^e (cos t - 3 sin t) - > 



* * -Bt * -Bt (26a 2 + ° l l ) 

x (t) = c e cos t + c„e sin t + r 

1 ( 6 + 1 ) 



(3a) 



in the region x 2 ^*-^ > c * an< ^ 



~ 3 1 — 3 1 

x (t) = -c e (B cos t + sin t) + c 0 e (cos t - 6 sin t) + a , 



2 

(3b) 



-St , , -St (26a 2 ta 1 ) 

x 0 (t) = c e cos t + c 9 e sin t - r 

1 ( 6 + 1 ) 



for X£(t) < c . Thus the trajectories in the phase plane are simply 
(arcs of) decaying (since B > 0) counterclockwise spirals, centered at 



(- 



(2Ba + a ) 

— — ) 

( 6 2 + 1 ) 



for the region x^ > c , and at 



( a 



2 



(2Ba + a ) 

— ) 

(6 + 1 ) 



(4) 



(5) 



3 



for the region < c . These centers will be asymtotically stable 
critical points if and only if they are actually above and below 

= c , respectively, e.g., (4) describes a critical point if and only if 

26 a 2 + a 

5 > c . 

(e z + 1) 

Note that, if and only if c = 0 , the origin is also a critical point, 
but obviously an unstable one* 

A representative view of the phase plane behavior can be obtained by 
considering the region x 2 > c . There, in the case 

(2Ba 2 + a ) < c(B 2 + 1) , 

the center for trajectory arcs is located in the region <_ c y with 

an x^ co-ordinate of -cl^ 0 • Thus, because the spiral arcs are 

counterclockwise, every arc emanating from a point in x 2 > c intersects 

the switching line at an x^ < 0 in some finite time, and exits the 

region. Furthermore, since the angular frequency is unity, then each 

arc spends no more time than tt in the region, with a time equal to tt 

occurring only when the center of the arcs is located on the switching 

2 

line (i.e., when (2£a 2 + a^) = c(8 4 * 1) ) . 

The situation when 

(2Ba 2 + a ) > c(6 2 + 1) 

is significantly different in that the spirals are now centered at a 
stable critical point actually located in the region x^ > c . Thus, 
some trajectories never exit the region, but decay down to the critical 
point, while others intersect the switching line and exit the region. 
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The "boundary" between the arcs that exit the region and those that decay 
to the critical point is the arc that terminates on the switching line, 



directly below the critical point, i.e., at 

( * 



( 6 ) 



This arc, which also exits the region, is denoted (*) in Figure 1. We 



let T denote the time represented by this arc, where clearly 



7T < T < 2tt . It is fairly easily shown that every other arc which exits 
o 



the region x 0 > c represents less time in the region than T * , i.e., 

/ o 

* 

T is the maximum time an arc can spend in x,> c and still exit there, 
o 2 

T q can be determined, using (x^ (t), x^ (t)) to denote the 
solution in the region x^ > c , by noting that the condition described 
by equation (6) and curve (*) of Figure 1 is eauivalent to 



* 1*0 ’ '“2 

x (T ) = x (0) = c . 
z o z 



(7) 



Using (3) and simplifying the resulting three equations, we arrive at: 

* 

-BT * * 

e = cos T -8 sin T , (8) 

o o 

which can easily be shown to have one and only one solution in (tt, 2tt) , 
* 

hence T^ is uniquely defined. 

The analysis of the region x^ < c is similar, and so is omitted. 

The only point of interest in this case is that, when, a critical point 

2 

is located in this region (i.e., + a^) < c(8 +1)) a "maximum 

time", T can be defined, similar to (7) by: 



x (T ) = a 
1 o 2 



x (T ) = x (0) = c . 
z o z 



(9) 
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S traighf orward calculation shows satisfies (8), hence, since the 

solution is unique, 



T = T * (10) 

o o 

(For the remainder of this discussion, we drop further reference to 
* 

T , replacing it by T .) Again we can show that no curve which 
o o 

exits x^ < c can spend more than time there. 
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Ill, THE CONDITIONS FOR A LIMIT CYCLE 

To derive the necessary and sufficient conditions for the existence 
of a limit cycle solution to (2), we begin with the solutions as given 
by (3). We let T^ denote the time spent by one trajectory arc in 
< c , and T^ the time on an arc in x^ > c . Then, since the 
system (2) is autonomous, and since trajectories must be continuous in 
the (x^,x^) plane, it follows that a necessary and sufficient condition 
for a limit cycle of period (T^ + T^) to exist is: 



X 1 (Tf) = x x (0) 



X 1 (°) = \ (t 2 ) 



x 2 (°) = x 2 ^ T i^ = c 

x *(0) = x *(T ) = c ^ 



T , T < T . 

1 2 — o 



( 11 ) 



where it is understood x^(0) represents x^(0 ) > etc. 

The condition given by (11) can be reduced by using the solutions 
given by (3) and simplifying the resulting expressions by eliminating 



c 2 and c 2 



This leads, in the case where neither T^ nor T 2 



equal tt , to the equations: 



cosh 6T- - cos T- . 

c i 1 + c 

1 sin T„ 1 



cosh BT 2 - cos T 2 
sin T 0 



= 0 



sinh BTj 
'1 sin T, 



- c. 



* sinh BT 2 
sin T„ 



= -2 



(1-6 )a 2 - 6a 1 

( 6 2 + 1 ) 



T l^ T o 
T 2 — T q 



( 12 ) 



(13) 



where 



* (2Ba ? + a 1 ) 

c - c = 2 Z 1 

C 1 1 



2 ’ 

(6 + 1 ) 



(14) 
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and the value of c for which this cycle occurs is given by 



c 



= 2 (C 1 + 




( 15 ) 



These equations express the necessary and sufficient conditions for 

existence of a limit cycle, except in the cases noted. (Note these later 

cases correspond respectively to c^ = 0 , if and only if, = tt and 

* 

c^ = 0 , if and only if, = tt .) In these cases, expressions similar 

•k 

to (12)- (14) can be derived. For example, in the case c^ = 0 , c^ ^ 0 , 
the condition for existence of a limit cycle of period (T^ + tt) is: 



sinh 3 tt 
cosh $7T + 1 



cosh gT 2 - cos 
sin T 2 



sinh 0T (1-0 )ot - 20a 

L + — t < x 

sin T 2 (20a 2 + a ) ’ 2 - o ' 



(16) 



We now develop an interesting, and powerful, geometric interpre- 

2 *2 

tation of the necessary and sufficient conditions when c^ + c^ >0 

Consider (12- (14), which represent the necessary and sufficient conditions 
* 

when c^c^ f 0 , and let 



F(t) 



sinh 0t 
sin t 



G(t) 



cosh 0 1 - cos t 
sin t 



(17) 



Then (12) - (14) is equivalent to the vector equation 
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G(Tp 




-g ( t 2 ) 


F(T 1 ) 


* 

" C 1 


F(T 2 ) 


1 




1 — 1 


M _ 




- — 



2 ( 26 a 2 + a^) 
(6 2 + 1 ) 



(18) 

0 1 

(1 - B 2 )a 2 - Ba 1 
(2Ba 2 + a x ) 

1 



if ( 23 + a ) ^ 0 . But this means that the three vectors in (18) are 

3 

linearly dependent, i.e., lie in the same plane in R . Clearly this 
plane cannot coincide with the z = 1 plane, yet since each vector in 
(18) terminates on the z = 1 plane, linear dependence then requires that 
the end points of the three vectors lie in a line in that plane. Thus, 
(12)-(14) can be satisfied if and only if the vectors 



G(T ) 




-G(T 2 )- 




0 
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, and 


(1 - 6 2 ) a - 6a- 


F(T 1 ) 




f(t 2 ) 




Z 1 

-(26a 2 + a 1 ) 



(19) 



lie in the same straight line. But by drawing the parametric curves: 

?! = { (G(x) , FCc)) 0 < x < T } , (20) 

1 — — o 

C 2 = { (-G (x) ,F(t)) 0 < x < T q } 

plus the point 

(1 - B 2 ) a - 6a 

P = (0, ) , (26a + a ) 4 0 , (21) 

0 (26a 2 + ap 

'fc 

we have that (provided ^ 0 ’ anc ^ (26a 2 + a^) 4 0) , any line 
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through P and both curves generates a limit cycle, and, conversely, 
o 

for any limit cycle, there must be a straight line through the point and 
both curves. 

Straightforward but laborious calculation, plus liberal application 

of hyperbolic identities, shows that both F(t) and G(*r) are positive 

and monotonically increasing on (0 , tt ) , and negative and convex down on 

(tt, 2tt) . Furthermore G(t) has its maximum to the right of 

(i.e. G(t) increases monotonically on (tTjT^)) , while F(t) has its 

maximum at a T < T . Also, and are convex up in the upper 

mo ± z 

half plane, intersecting at (0,3 ) , and convex down in the lower half 
plane, being asymptotic to the lines 



and 



sinh 3^ 0 

U cosh 3 tt + 1 w 



( 22 ) 



.sinh Btt _ 

cosh 3 tt + 1 W 



(23) 



respectively, as t tt (w-*-±°°) . Thus and must have 
the approximate shapes shown in Figure 2. Figure 3 shows one situation 
which would lead to a limit cycle. Observe from (12) that the ratio of 
the horizontal intersections on and > which we denote: 

* 



G(T X ) 

g(t 2 ) 



(24) 



satisfies (from (14) and (15)); 



c 



1 

2 



(c l + 



(23a 2 + 

(3 2 +1) 



(1 + P) 



(1 ~ p) 



(25) 



10 



Thus, since p is just the ratio of horizontal co-ordinates, we could 
easily compute, from the geometrical configuration, the value of c for 
the cycle shown in Figure 3. (Assuming values for B^ , and 

are known.) 

Equations (24)— (25) will be especially useful in discussing the 

range of values of c for which cycles can occur, and, since it will 

become important later, note that on any interval not including p= 1, 

c as given by (25) is a monotonically increasing [decreasing] function 

of p for ( 2 go <2 + a^)> 0 [(230^ + a )< 0] 

In the discussion so far, this geometrical interpretation has been 

limited to the cases c^c^ ^ 0 > and ^Ba^ + a^) ^ 0 . However, the 

2 *2 

interpretation can easily be extended to all cases where c^ + c^ >0 
by extending the interpretation of a line intersecting the point and both 
curves to allow intersection in the limit at infinity. Thus, for example, 
the case c = 0 , c^ ^ 0 and (230^ + a^) ^ 0 is treated (refer to 
(24)) as the limiting case 

G(T ) + ± » . 

Now, referring to Figures 2-3, it is obvious “intersection” in this 
limiting case occurs if and only if the straight line is through P q , 
parallel to the asymptote of £ , intersecting ^ • (Figure 4). That 

is, there must be a T^ such that 



(1 - B 2 ) a - 6a 

F(T 2> + ■ ( 26 ^ +a - 2 r 

- G(T 2 ) 



sinh Btt 

cosh Btt + 1 



(26) 



But referring to the definitions of F(t) and G(t) it can be seen (26) 



11 



is identical to (16). A similar process can be followed when 
c 1 * = 0 , c^ ^ 0 , and (23a 2 + a ± 0 . 

Lastly, interpreting (19)-(21) where (23a^ + a^) =0 as the 
limiting case when P q approaches the point at infinity along the verti- 
cal axis, we see that an intersection" can occur if and only if the 
points on and ^ li e on the same vertical line, or equivalently 

G(T^) = - G(T 2 ) • However, since by (14)- (15) this case corresponds to 

* 

c l= c ± = c , this interpretation immediately satisfies (12), and yields 
from (13) a cycle occuring for every value of c satisfying: 



c 




sinh 3T^ 
sin 



sinh 3T 2 
sin T 2 



(27) 



where and are solutions of G(T^) = ~G(T^) , and the absolute 

value of the quantity in brackets represents the vertical distance 
between the two points of intersection. 

Thus, including the limiting cases, we see that for all 

2 

c l + c i > 0 , the geometrical interpretation is equivalent to the 
solution of (12)- (14), and hence it expresses the necessary and sufficient 
condition for existence of a limit cycle. 
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IV. EXISTENCE OF LIMIT CYCLES 



In this section we apply the geometric construction developed above 
to derive necessary and sufficient conditions under which a limit cycle(s) 
will exist. Our main conclusion will be to determine existence conditions 
in terms of 3 > and • Specifically, we show a bound, depen- 

dent solely on 3 , such that no cycles exist for any values of c if 
the ratio [a^/a^] exceeds this bound, and, for all values less than this 
bound, there is a c > 0 , depending on 3 , a. , and such that 

a cycle will exist if and only if c < c 

J J ii— max 

We begin by noting the following result, alluded to in the introduc- 
tion to this paper: 



Lemma 1: A limit cycle will exist for values of c > 0 if and 
only if a cycle of the same period exists for - c 



The proof of this, which follows straightforwardly from the oddness 
of all functions in (2), is omitted here. 

The following theorem gives the condition under which no cycles of 
any period are possible, irrespective of the biasing. 



Theorem 1: 



There is an n (3) such that (2) will have no 
max 



limit cycle solutions if either 



(a) c *2 = 0 , or 



(b) 



and 



> n (8) 

max 



Proof: Let M(B) = - max 



(tt> 2tt) 



sinhBt 



sin t 



> 0 , 



(28) 
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and 



n 

max 



( 8 ) 



1 - 8 - 28M 

(6 + M) 



> -28 



(29) 



(Note - M is the maximum vertical coordinate on both £ and 5 in 
the lower half plane.) It is then easily shown that conditions (a) and 
(b) are equivalent to 



- M < 



(1 - 6 2 )a 2 - 30^ 

( 28 a 2 + a 1 > — ^ 



However, this implies that P q (as given by (21)) lies above the line 
joining the maxima of and ^ in the lower half plane, but at or 

below the juncture of and ^ in the upper half plane (Figure 5). 

The convexity of and 5^ and the quandrants in which the curves 

are located immediately imply no straight line through P q can intersect 
both ^ and ^ > hence no cycle can exist. | 

Note that, by defining 

, a 2 4 0 , (30) 

P q can be equivalently defined 




P Q (B,n) 




(1 - 8 2 ) - 8ti 
(23 + n) 



(31) 



Then, if the situation when = 0 is considered as the limiting case 

n "*■ ± 00 , the previous theorem implies cycles can exist only if 

— oo < n < n (3) • We now complete the result by showing that if 
— max 

” 00 < n £ T) (3) » a limit cycle will exist for at least one value of c 



14 



Theorem 2: If n = n (B) , a cycle will exist if and only if 

max J 

c = 0 . 

Proof: If n = n (6) , then 

max 

P q = (0,-M) , 

that is P lies on the line joining the maxima C- and £ in 
o 12 

the lower half plane, (Figure 5). But the convexity of and 

imply this is also the only line through P intersecting both 
2 o 

curves. Thus one and only one cycle will occur in this case. From 
Lemma 1, c = 0 . | 



The above result is the reason why, in the introduction, we commented 

that existence of a limit cycle in the unbiased case was a necessary and 

almost sufficient to give cycles for some biasings, for, as the following 

shows, only when n (6) = n can a cycle exist only for c = 0 

max 

Theorem 3: If - 00 < n < n (6) , then there is a c > 0 , 

max max 

such that a limit cycle will exist if and only if 



Proof : By Lemma 1 we restrict ourselves to c >_ 0 . The actual 

proof is divided into three subcases; -26 < n < n (6) , 

max 

n = - 26 , and - 00 < n < -26 

Case I: The case -26 < n < n (6) > corresponds (with our re- 

max 

striction a ^ 0) to (260^ + ot^) > 0 , and so, by (25), c > 0 

corresponds to -1 <_ p < 1 . Also, from (29) - (31), we can show 

that in this case P (n>B) is located in the lower half plane, 

o 

below the line joining the maxima there of anc ^ * However, 
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TT 

with P so located, we can define 0 , 0 < 0 < — , as the 

° ° ° 2 
largest clockwise angle from the horizontal for which a straight line 

through P intersects in the third quandrant. Then clearly 

o i 

any line through P q with slope, , will intersect both an ^ 

C 2 if and only if 

| m | <_ tan 0^ . (Figure 6) 

However, since p must vary continuously with m^ (along lines 

intersecting both curves), and since p can approach 1 only if the 

intersecting line approaches the vertical, then p must be bounded 

away from 1 , and there exists a p such that intersection will 

J max 

occur for c > 0 if and only if - 1 < p < p < 1 . It follows 

— max 

from (25) that a cycle will exist for c >_ 0 if and only if 



(2ga ? + a ) 

0 < c < c » — — 

- " (g 2 + x) 



(1 + p 

(1 - P 



) 

max 

max^ 



> 



since the right hand side of (25) is strictly increasing for 

(2Ba„ + a, ) >0 and -1 < p < 1 . Clearly c depends on 3 , a. 

2 1 — max l 

and <*2 

Case II : When -23 = n , (23ot 2 + a^) = 0 , and as we have discussed 
before in this case limit cycles exist if and only if c satisfies 
(27), where in (27), the quantity: 

sinh 3T sinh 3T 2 



represents the distance between £ and £ on any vertical line 
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intersecting both (since P q C n ) is imagined passed to a limit at 

infinity along the vertical axis). Since the distance must clearly 

have a non-zero minimum, and since (27) involves its reciprocal, 

it follows straightforwardly that there is a c such that c , 

max 

as defined by (27) satisfies 

0 < c < c 
— — max 



Case III : When - °° < n < -23 , the proof is essentially similar 

to case I, except P (h,3) is located in the upper half-plane, 

o 

above the juncture of and ^ • | 



Theorem 3 completes the description of the necessary and sufficient 
conditions under which a limit cycle (s) will exist. Note that although 
this theorem predicts at least one cycle for every value of c satisfy- 
ing 0 < c < c , it does not address the possibility of more than 
— 1 1 — max 

one cycle occuring at a single value of c . We shall pursue this in 

the next section. However, we first investigate in some more depth the 

properties of n (3) 
max 

Elementary methods show that, for a given 3 > M(3) as defined by 

(29) is given by 



where T (3) 
m 



M(3) = - ‘ 

is the solution, in 



sinh 3T (3) 

m 

sin T (3) 
m 

(tt,2tt) , of the equation 



tan t = “ tanh 3t 

p 

Since the right hand side of this equation can be shown to be monotoni- 
cally decreasing in 3 (for fixed t ) , it follows that T m (6) is 
monotonically decreasing in 3 > and 
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lim T (3) = tt 

n V m 



Further direct calculation will show 



, m/oN T (6) cosh 3T (3) 
dM(3) _ _m m 

d3 ~ sin T (3) 



> 0 



m 



since tt < T . Thus M(B) is a monotonically increasing, and clearly 
m 

M(B) 00 as $ ^ 00 

With M(B) thus fully characterized, it is easily shown that 

ri (B) » as given by (29), is a monotonically decreasing function of 
max 

B , asymptotic to the curve r| = -2B as B 00 . Figure 7 shows the 

general dependence of r > max > and it is fairly obvious that 

n (3) = -23 , 3 > 1 . 

max 
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V. BIFURCATION OF THE LIMIT CYCLES 



The previous section showed that for - 00 < H Ji max(3) , there was 

a c (n,3) such that at least one limit cycle would exist for every 
max J 

c satisfying 0 < |c < c . In this section we complete the 

description of the limit cycles by showing there are values of c that 
allow more than one cycle, how these values are determined, and how the 
transition from single to multiple cycles occurs. 

Consider the situation represented by figure 8. When the straight 
line intersects two curves and ^ > the horizontal co-ordinates of 

the intersections and slopes there of these curves are denoted 
W 1 * m i anC * w 2 » ra 2 respectively, and a is the angle between the 
line and the horizontal. Then, using elementary calculus, it can be 
shown that: 



{ — > = — { 
da w ^ w ^ 



m 



+ 1 } 



(m 2 -m 1 ) 



(m -m ) (m -m ) 
la z a 



However, referring to (24), this equation becomes 



it, . 2 . tyV 

da p m a ^ (m -m ) (m -m ) 

la z a 



(33) 



where m^ is the slope of the intersecting line. This equation allows 
us to determine, relatively easily, the points where changes sign. 

But from (25), c is a monotonically increasing function of p for 
(23 + n) < 0 , and a monotonically decreasing function for (23 + n) < 0 
Thus every time changes sign (for fixed values of 3 ,n ) , a bifur- 

cation occurs, i.e. two cycles exist for the same values of c 

With this formulation, and our main intent being to characterize 
the algebraic sign of as given by (33), we proceed. As will emerge 
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more clearly later in the discussion, there are several special subcases, 

of increasing complexity, that must be considered. We shall not do this 

exhaustively, but essentially look in full detail only at one of these 

cases, then comment on the others. In essence, the complex cases arise 

when P is located "near” (0,-M) , so that a single line through P 
o o 

can intersect (or twice in the same quadrant. 

Consider the case where P q is located in the lower half plane, 

and at or below the intersection of the tangent to at T with the 

1 o 

verticle axis. (Figure 9). Since we can show m . (T ) = -1 , this condi- 

1 o 

tion can be reduced, using the definition of T , to 



where 
tan 0 



0 < 20 + n £ - 

F(t) was defined in (17). Note 
> 1 . 



(1 + B 2 ) 

2 (F (T q ) " B) 

also this immediately implies 



o 

In this situation, due to the symmetry and convexity of and 

, no line thru P can intersect either curve more than once in any 
2 o 

quadrant. Furthermore, as noted in the previous section, since 

(20 + n) > 0 it follows cycles with c >_ 0 occur only for 

-1 £ p < 1 , which equates to lines thru P q whose slopes satisfy 

- tan 0 < m < 0 , 

O — OF — 

where 0 was defined in Theorem 3 and Figure 6 . 
o 

The easiest way to arrive at the information needed to determine 

is by describing the variation of m and m^ as a increases from 

-0 to zero. In this case, observe that 0 is determined by the line 
o o 

thru P to (G(T ), F(T )) , the end point of . Thus as a 

o o o 1 

increases to zero, m^ increases monotonically (due to the convexity of 
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from -1 to a positive value less than 

sinh Btt .. 

cosh 3 tt + 1 

Since, as noted above, tan >_ 1 , then as a increases to zero, 

m = tan a must increase monotonically from -tan to zero. Note 

that since a line through P as located here cannot intersect 

o *1 

tangentially (except perhaps at the endpoint) , it follows 

m < m_ for -0 < a < 0 

a 1 o 

Las t ly , since m < -1 at -0 , then the line through P at -0 

a — o o o 

must intersect £ ^ second quadrant. (Since here asymtotic 

to a curve whose slope is greater than -1). Thus as a increases, m^ 

will initially monotonically decrease from a value in (-1, 0) to 

sinh 3 tt 
cosh 3 tt + 1 * 

at which time the intersection point transitions to the fourth quadrant 

(and also p passes through zero, changing sign), then m^ increases 

monotonically to a value at a = 0 which is the negative of m^ at 

a = 0 . It thus follows m^ < m^ < 0 in second quadrant, and 

m^ < < 0 in the fourth quadrant for a <_ 0 

If m^ , m^ and m^ are now plotted based on the above analysis, 

a picture like Figure 10 will emerge. Observe that and must 

intersect exactly once for -0 < a < 0 . Now if we denote the inter- 

J o 

section of m, and m^ as -a, , and that of m, and m as -a , 

12b 1 a t 

where necessarily -a, < -a , we can easily construct Table 1 to give 

b t 

the algebraic signs necessary to interpret (33). 
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TABLE I 



-0 < a < - a, 

o b 



-a, < a < -a 

D t 



-a < a < 0 



m l ■ 


m 

a 


+ 


+ 


+ 


m 2 ' 


m 

a 


+ 


+ 


- 


m 2 




+ 


- 


- 


P 




+ 


+ 


- 


d£_ 

da 




+ 


- 


- 



From this table it is clear that multiple cycles (for the same value 

of c ) can exist in some cases. Furthermore this multiplicity can be 

of order exactly two, since only changes sign once. (Again recall 

c is a monotonic function of p .) Lastly, the multiple cycles can 

occur only in a band of values of c whose lowest bound is strictly 

positive. In fact, due to the sign of ^ on -0 < a < -a, , the value 

da o b 

of c at which the multiple cycles first occur is precisely the one 

found from the value of p for the line with slope -0 . It should be 

o 

noted that since the phase plane is two dimensional, when double cycles 

occur at least one will be unstable. This behavior, existence of a 

single cycle for 0 < c < c, , then two cycles for c, < c < c , seems 

b b — max 

best described by the statement that a bifurcation occurs at c^ 

The situation when (26 + ri) £ 0 can be easily constructed in a 
manner to the above, and yields essentially the same result - a unique 
cycle will exist for 0 < I cl <c, where c, < c , there for 
c^ <_ | c | <_ c exactly two cycles occur. 

The other possible situations which arise become progressively more 
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complicated from a bookkeeping point of view. The reason for this lies 

in the geometry, with the convex nature of the curves involved. For, 

since is convex, the curve formed by plus the extension of the 

tangent line at the endpoint (at T ) to the intersection with the 

o 

verticle axis is also convex. Thus, where 

. ik+jh. < 26 + n , 

2 (F (T ) -3) 
o 

P lies "outside" this extended curve, and hence for at least some 
o 

angles o , a line through P q with slope m^ will intersect twice. 

Thus m^ cannot be described by a single-valued function of a , and 

the analysis becomes more complex. However, as long as P^ is located 

below the line joining the endpoints of an d ^ t ' ie third and 

fourth quadrants, the same basic behavior as was found in the simple 

cases above can be shown to exist. 

However, when P lies on the line joining the endpoints of C-, 
o 1 

and C , i.e. when 

1 2 

o < 26 + 1 * - pVt Z 6 . (34) 

o 

a qualitative change in the behavior occurs. For at this point, as 

Figure 11 shows, there are two solutions (points (l)-(3) and (2)-(4)) for 

which p = -1 , and hence c = 0 . Thus, in addition to the bifurcation 

at c^ > 0 , there is now a second bifurcation at c = 0 . As P q 

continues to move toward (0,-M) , eventually a condition is reached 

where multiple cycles exist at all values of c 

Thus, in conclusions, we have shown that, in all cases, there is a 

c such that two cycles exist for c, < I cl < c , and in certain cases 
b b — 1 1 — max 
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there is also a 



0 < 




1 

C b 



- °b 



such that two cycles can exist for 
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VI . SUMMARY 



In this paper we have considered limit cycle solutions to 
x 1 = -23x 1 - (3 2 +1)x 2 + sgn (x 2 -c) 

x 2 = X^^ + a 2 sgn (x 2 ~c) 

We have derived the necessary and sufficient conditions for limit cycle 
solutions to exist, and have presented the conditions in terms of an 
equivalent, easily visualized, geometric condition* With the use of this 
geometric interpretation, we then showed that, with one exception, 
existence of a limit cycle for c = 0 is both necessary and sufficient 
for cycles to exist for some range of c 4 0 . Lastly, we have shown 

that multiple cycles, for some values of c , are possible is essentially 
all cases* 
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Appendix 1 



We consider the system 



v = A v + sgn (rv^ + sv^ - c) 



(A-l) 



where u is a constant vector, r + s >0, and A is a real, 

— o = 

2 

constant 2x2 matrix whose characteristic polynomial is X + 2yX + 
2 2 

(y 4- to ) , where y, a) > 0. First observe that the non-singular 
transformation 



S ■ B v, 





reduces (A-l) to 



-1 



w = (BAB )w + B u sgn 



(a)W 2 - c) 



(A-2) 



where, because of similarily, BAB has the same characteristic 
polynomial as A, i.e. 



x 2 + 



2y\ + (y + 



o> 2 ) 



Since to > 0, both off-diagonal terms in BAB’*' must be non-zero, 
and so BAB'*' will have the form, 
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BAB 



-1 



11 



21 



12 



22 



where d 2 ^d 12 * °* Thus 4 0, and hence the transformation 



X = c w , 
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— 




is non-singular. This transformation reduces (A-2) to 



± = (CBAB 1 C 1 )v 



+ C B _u^ sgn (uy 2 - c) 



(A-3) 



where 



C (BAB 1 ) C _1 



r 



( d 2 1 + d 22 ) 



^ d ll d 22 d 21 d 12' ) 



Again by similarity, this matrix has the characteristic polynomial 

2 2 2 
A + 2y\ + (y + u ) . 



By computing the characteristic polynomial directly, we see 



-2y = (d 



11 



d 22 ) 



2 

Y 



+ 03 



d ll d 22 d 21 d 12 
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Thus , 



2 2 

y x = - 2 X y x - (y + <*> )y 2 + ^ sgn (<oy 2 - c) 



(A-4) 



y 2 = y l 



+ u sgn (toy - c) 



where 



u = C B u 



But, now the time scaling oat = t reduces (A-4) to 



y 1 = ~2(^)y 1 - ((^-) 2 + 1) (coy 2 ) + — • sgn ((toy 2 ) - c) 



(wy 2 ) = y x 



+ u 2 sgn ((toy 2 ) - c) 



which, with the transformation 



x., = 



y l • x 2 ' " y 2 ' S ‘ 1 



a. = — , a„ = u 
1 to 2 2 



becomes ( 2 ) . 
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